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1.0  Introduction  and  Overview 


The  purpose  of  this  effort  has  been  to  provide  ITT  Systems  &  Sciences  with  a  knowledge 
based  tracking  capability  that  will  support  a  space  time  adaptive  processing  (STAP) 
environment.  This  has  included  the  development  of  a  basic  multiple  target  tracker  and  the  design 
and  testing  of  several  knowledge  based  rules.  A  rule  book  containing  25  potential  knowledge 
based  rules  was  developed  and  is  presented  in  Volume  V. 

For  the  purpose  of  ITT’s  application,  the  main  elements  of  the  tracker  software  can  be 
imbedded  in  a  larger  STAP  simulation  by  removing  the  GUI.  The  key  tracking  modules  are 
setup_Tracker_6  and  Run_Tracker_8. 

The  Run_Tracker_module  allows  the  use  of  three  types  of  tracking  filters.  These  include: 

1.  an  uncoupled  two  state  alpha  beta  filter  with  position  and  velocity  component 
states, 

2.  an  uncoupled  three  state  Kalman  filter  with  position,  velocity,  and  acceleration 
component  states,  and 

3.  an  extended  four  Kalman  filter  with  both  x  and  y  position  and  velocity  component 
states. 

The  first  two  filters  use  a  one  dimensional  measurement  vector  containing  the  report 
position  component.  The  extended  Kalman  filter  uses  a  three  dimensional  measurement  vector 
containing  x  and  y  report  position  and  “pseudo”  Doppler,  the  latter  defined  as  range  times  range 
rate. 


The  tracking  software  as  delivered  has  been  imbedded  in  a  GUI  structure  that  makes  it 
easy  to  exercise  the  tracker  under  a  variety  of  conditions.  Using  the  GUI,  the  user  can 
interactively  select  existing  or  new  target  scenarios  and  tracking  options  prior  to  tracker 
execution.  Figure  1-1  shows  the  list  of  available  options. 
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GUI  Interface 


•  Default 

•  Existing 

•  New 


•  Radio  frequency 

•  Range  resolution 

•  Azimuth  resolution 

•  Doppler  resolution 

•  CPI’s  per  dwell 

•  False  Alarm  Probability 

•  Target  SNR 

•  Pulse  Repetition  Frequency 


•  Default 

•  Track  Parameters 

•  Promotion  Logic 


•  LP/HP*  Hit  Table 

•  LP/HP  Miss  Table 


*  LP  and  HP  denote  low  priority 
and  high  priority,  respectively 


•  Track  Filter  Type 

•  LP/HP  Gate  Multiplier 

•  LP/HP  Tentative  Gate  Radius 

•  Maneuver  Noise 

•  Alpha/Beta  Filter  Gains 

•  Target  Priority 
(max  of  5  targets) 

•  Target  Turn  Interval 
(max  of  5  targets) 

•  Tracker  Gate  Type 


Figure  1-1:  Tracker  GUI  Options 

2.0  Tracker  Description 

The  first  part  of  this  section  contains  a  narrative  description  of  the  basic  multiple  target 
tracker.  Detailed  documentation  has  been  provided  in  Section  3.  This  is  followed  by  a 
discussion  of  three  knowledge  based  tracking  rules  that  can  be  used  to  support  the  STAP 
processor.  Simple  tracking  simulations  have  been  provided  to  illustrate  the  use  of  these  rules. 

2.1  Basic  Tracker 

A  high  level  description  of  the  multiple  target  tracker  function  is  shown  below  in  Figure 
2-1.  The  tracker  processes  reports  on  a  per  scan  basis  and  makes  use  of  a  track  table  that 
contains  several  track  attributes  as  well  as  kinematic  information.  The  key  attributes  include  the 
tracks  identification  number,  its  state  value,  which  is  a  measure  of  track  quality,  and  its  status. 
Track  status  can  be  dropped,  tentative,  or  firm.  In  the  current  software,  a  dropped  status  is 
assigned  to  a  track  whose  state  has  been  demoted  to  zero,  a  tentative  status  is  assigned  to  a  new 
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track  formed  by  an  unused  report  with  a  state  value  initialized  to  one,  and  a  firm  status  is  given 
to  any  track  with  a  state  value  greater  than  one.  The  tracker’s  function  consists  of  a  correlation 
section  together  with  association  and  track  maintenance  sections. 


Figure  2-1:  High  Level  Tracker  Flow 

The  correlation  section  is  performed  by  testing  for  inclusion  each  current  report  against 
each  extrapolated  track  gate  and  forming  a  binary  correlation  matrix  with  ones  in  the  capture 
positions.  Tentative  track  gates  formed  from  unassociated  reports  are  extrapolated  by  centering  a 
circle  of  kinematically-determined  maximum  radium  about  the  report,  whereas  firm  track  gates 
are  extrapolated  by  using  the  track  filter  prediction  equations.  Track  gate  sizes  are  typically  a 
function  of  both  the  measurement  and  prediction  uncertainties,  as  well  as  the  track’s  maneuver 
status.  A  distance  matrix  of  the  same  size  as  the  correlation  matrix  is  also  formed  containing  the 
report-to-gate  center  distances. 

The  association  and  maintenance  section  uses  the  above  correlation  and  distance  error 
arrays  to  assign  reports  and  tracks  in  a  unique  manner.  In  the  even  that  multiple  reports  are 
common  to  a  given  track  gate,  or  multiple  tracks  capture  a  common  report,  the  association  logic 
will  resolve  the  conflict.  Furthermore,  any  reports  Aat  are  not  assigned  to  an  existing  track  will 
be  used  to  spawn  a  new  track  designated  to  have  a  tentative  status.  Two  simple  association 
logics  are  currently  available.  There  is  a  “nearest  neighbor”  logic  that  assigns  the  closest 
captured  report  to  a  given  track  and  there  is  a  “venerable  track”  logic  that  associates  the  oldest 
common  track  to  a  given  report. 
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After  all  tracks  and  reports  have  been  associated,  the  track  state  promotion  logic  is 
applied.  Track  states  are  updated  using  either  a  “hit”  or  a  “miss”  table,  depending  on  whether 
they  correlate  with  a  report.  These  state  tables  are  usually  designed  to  require  specified  numbers 
of  successive  hits  and  misses  before  a  track  is  declared  as  firm  or  dropped.  Two  examples  are 
given  below.  In  Case  1,  tables  Hitl  and  Missl  allow  a  tentative  track  to  build  up  to  firm  status  at 
state  level  4  after  three  successive  this,  and  to  demote  with  each  miss  until  it  reaches  state  0 
where  it  is  dropped.  In  Case  2,  tables  Hit2  and  Miss2  show  a  more  complex  strategy.  A 
tentative  track  promotes  up  to  a  firm  status  of  3  after  two  hits.  However,  there  are  hold  states  4,5 
and  6  which  allow  the  track  to  recover  its  firm  status  more  quickly  after  a  single  miss. 


Case  1 

State 

1 

2 

3 

4 

Hitl 

2 

3 

4 

4 

Missl 

0 

1 

2 

3 

Case  2 

State 

1 

2 

3 

4 

5 

6 

Hit2 

2 

3 

3 

5 

3 

3 

Miss2 

4 

6 

5 

0 

4 

0 

Once  track  state  logic  has  been  applied,  the  track  table  is  updated.  All  tracks  with  state 
values  of  0  are  removed  from  the  table.  Furthermore,  all  reports  that  were  not  associated  with 
extant  tracks  are  used  to  spawn  new  tentative  tracks  that  are  added  to  the  table. 

2.2  Knowledge  Based  Rules 

Knowledge  based  rules  make  use  of  extended  map  and  intelligence  information,  and  are 
used  to  improve  the  tracker’s  ability  to  support  STAP  processing  requirements.  Key  issues 
include  the  tracking  of  targets  in  regions  adjacent  to  large  discretes,  and  in  shadow  zones  that  are 
blocked  from  the  radar-s  line  of  sight.  Using  map  information  the  tracker  should  also  be  able  to 
anticipate  a  target  maneuver  that  will  be  required  to  avoid  obstacles.  Queued  with  this 
information,  the  tracker  will  apply  appropriately  shaped  track  gates  that  enhance  its  capability  of 
capturing  reports  while  maintaining  reasonable  gate  sizes.  Four  rules  are  discussed  below,  along 
with  simple  tracking  simulations. 

In  the  following  examples  single  target  tracks  are  displayed  graphically  with  the 
following  conventions: 
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2.2.1  Tracking  Legend: 

+  Predicted  gate  center  position 

*  Measured  report  position 

o  Coasted  track  position 

•  Extrapolated  position  of  a  missed  detection 

d  Dropped  tentative  track 

D#  Dropped  firm  track 

#  Corresponds  to  age  (scans)  of  dropped  track 

In  addition,  all  the  results  were  obtained  using  an  uncoupled  Kalman  tracking  filter  and 
all  simulations  assumed  a  ten  second  scan  period. 

2.2.2  Maneuver/Obstacle  Rules: 

Both  alpha  beta  and  Kalman  tracking  filters  do  a  good  job  with  targets  that  move  along  a 
constant  heading  with  a  fixed  speed.  Deviations  from  a  straight  path  cause  prediction  errors  to 
occur  and  can  ultimately  result  in  a  dropped  track.  Therefore  it  is  important,  whenever  possible, 
to  anticipate  target  maneuvers  by  several  scans.  This  allows  time  to  make  such  adjustments  as 
increasing  the  gate  size  and  track  gain,  or  using  a  shaped  gate  to  allow  for  across  track  deviations 
caused  by  target  turning. 

If  a  track  approaches  an  obstacle  whose  across  track  extent  is  H,  a  maneuver  can  be 
anticipated  to  occur  within  a  time  extent  no  longer  than  Tmax-  For  this  discussion,  assuming  a 
constant  target  speed  v  and  a  maximum  possible  acceleration  Amax,  this  extent  is  given  by: 

V 

where 


denotes  the  radius  of  curvature  of  the  target  turn  required  to  clear  the  obstacle.  Somewhere 
within  this  time  period  the  tracker  should  apply  its  maneuver  logic. 

Figures  2-2  to  2-4  use  the  same  section  of  track  to  illustrate  the  effect  of  different  gating 
strategies  on  the  tracker’s  ability  to  handle  a  maneuvering  target.  Performance  is  computed  for  a 
constant  speed  (250  meters/sec),  high  SNR  (20  dB  at  mid  range)  Swerling  1  target,  as  it 
approaches  an  obstacle  (shaded  rectangular  region)  from  below  and  begins  to  turn  after  the 
twenty  third  scan. 

Figure  2-2  shows  the  response  of  an  uncompensated  tracking  filter  using  a  standard  track 
gate  centered  on  the  predicted  gate  center,  and  oriented  along  and  across  range  with  a  size 
dependent  on  both  the  measurement  and  prediction  errors.  No  target  acceleration  has  been 
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assumed  and  no  maneuver  noise  has  been  added.  While  the  straight  line  section  is  handled 
adequately,  the  initial  track  begins  to  lose  the  target  after  the  turn  begins,  whereupon  it  is 
demoted  to  a  dropped  status  on  the  next  four  scans.  The  three  kilometer  circular  gates  indicate 
newly  spawned  tentative  tracks  that  were  subsequently  updated  to  firm  status. 
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Figure  2-2:  Uncompensated  Tracker  Response  for  Maneuvering  Target 
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The  first  maneuvering  target  rule  specifies  the  use  of  shaped  elliptical  gates.  Figure  2-3 
shows  the  same  target  as  tracked  using  a  combination  of  two  gates,  both  centered  on  the 
predicted  gate  center.  One  gate  is  oriented  along  and  across  range,  with  a  size  determined  by 
range  and  angle  measurement  uncertainties.  The  second  gate  is  oriented  along  and  across  the 
targets  instantaneous  track  and  its  size  is  a  function  of  the  maximum  turning  acceleration  of  the 
target,  assumed  here  to  be  0.5g  units.  Reports  capture  occurs  if  the  measurement  falls  within 
either  gate.  Except  for  a  missed  report  early  in  the  linear  part  of  the  trajectory,  causing  the  gate 
to  swell  initially  and  then  settle  down,  the  track  is  maintained  throughout  the  maneuver.  Each 
lower  case  d  symbol  indicates  a  dropped  tentative  track.  This  occurs  when  such  a  track  fails  to 
capture  a  report  on  the  following  scan.  The  large  three  kilometer  circles  denote  tentative  tracks 
that  were  successfully  promoted  to  firm  status.  Finally,  the  dots  occurring  in  both  Figures  2-2 
and  2-3,  shown  extrapolated  from  the  straight  line  section  of  the  approaching  track,  represent 
missed  detections.  While  actual  gates  existed  for  these  cases,  they  have  been  drawn  here  only  or 
situations  in  which  reports  were  captured. 
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Figure  2-3;  Compensated  Tracker  Response  for  Maneuvering  Target. 

Track  Oriented  Elliptical  Gates  used  with  0.5g  Across-Track  Acceleration  Allowance 


t 
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The  second  maneuvering  target  rule  specifies  the  use  of  a  gate  whose  shape  is  determined 
solely  by  centripetal  turning  mechanics.  Figure  2-4  shows  the  target  being  tracked  using  a 
“smile  shaped”  gate  whose  shape  is  determined  by  the  kinematic  constraints  imposed  on  a 
constant  speed  target  undergoing  a  centripetal  maneuver.  Let  t  range  over  the  time  interval  from 
the  last  track  update  till  the  next  predicted  report  acquisition,  Tscan  seconds  later.  The  maneuver 
envelope  is  given  by  the  xy  locus  of  points,  oriented  along  and  across  the  track,  and  generated  by 
the  equations: 


x  =  vt  +  psinO 
y  =  p(l-cos^) 

0  =  v(T^,„,-t)l  p 


where  the  radius  of  curvature  p  is  defined  as  above.  For  this  gate,  a  track  capture  occurs  if  the 
measurement  ellipse,  centered  on  the  measured  report  and  oriented  along  and  across  range, 
intersects  the  smile  locus.  As  in  the  previous  example,  this  gate  is  able  to  maintain  the  track 
throughout  the  target  maneuver.  One  advantage  of  this  gate  is  that  is  has  a  relatively  small  area 
as  compared  with  other  maneuver  gates  and  this  makes  it  less  likely  to  capture  any  false  alarms. 


x-coordinate  (km) 


Figure  2-4:  Compensated  Tracker  Response  for  Maneuvering  Target.  “Smile  Shaped 
Centripetal  Gates  used  with  (0.0  -  0.5)  g  Range  of  Across-Track  Acceleration  Allowance 


8 


2.2.3  Shadow  Rule: 


The  shadow  rule  provides  a  means  of  preserving  firm  tracks  that  enter  regions  shadowed 
from  the  radar  line  of  sight.  If  the  predicated  track  gate  center  falls  within  a  designed  shadow 
region,  both  the  track  state  and  gate  size  will  be  frozen.  Upon  emerging  from  the  shadow,  state 
promotion  resumes  and  the  gate  size  will  not  be  allowed  to  exceed  a  maximum  value.  In  the 
tracker  scenario  used  in  Figures  2-5  and  2-6,  only  a  minimal  amount  of  acceleration  noise,  0.05  g 
units,  was  used  in  order  to  maintain  a  straight  line  coasting  of  the  track  through  the  shadow.  As 
previously,  the  target  has  a  speed  of  250  meters  per  second  and  the  update  scan  time  is  10 
seconds. 

Figure  2-5  shows  a  section  of  tracker  response,  when  no  shadow  rule  is  in  effect,  for  a 
target  approaching  a  shadow  zone  (shaded  rectangle)  from  below.  The  three  dots  denote 
extrapolated  positions  of  the  initial  track  where  no  reports  were  captured.  After  four  demotions 
the  firm  track  that  initially  entered  the  shadow  was  dropped,  as  indicated  by  the  D14  symbol.  A 
new  tentative  track  was  spawned,  and  promoted,  once  the  predicted  gate  positions  moved  beyond 
the  shadow.  Note  that  the  D  symbol  has  been  placed  at  the  last  updated  track  position,  just  prior 
to  entering  the  shadow,  where  the  track  was  14  scans  old. 


Figure  2-5:  Tracker  Response  for  Target  Rying  through  Shadow 
No  Shadow  Rule  Applied 


Figure  2-6  shows  the  track  history  for  the  same  target-shadow  scenario  when  the  shadow 
rule  was  in  place.  The  circles  in  the  shadow  denote  coasted  track  positions  at  which  the  track 
state  was  held  fixed.  Once  the  predicted  gate  position  emerged  from  the  shadow,  there  was  a 
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moderate  increase  in  gate  size,  after  which  it  settled  down,  and  the  original  firm  track  continued 
undisturbed. 


Figure  2-6:  Tracker  Response  for  Target  Flying  through  Shadow 
Shadow  Rule  Applied 


2.2.4  Discrete  Rule: 

By  tagging  large  radar  returns,  or  discretes,  the  STAR  processor  can  exclude  regions 
containing  them  from  its  covariance  matrix  element  formation  and  thereby  not  use  up  limited 
degrees  of  freedom  on  their  cancellation.  The  discrete  rule  allows  the  tracker  to  coast  through 
any  region  containing  one  of  these  tagged  returns  and  to  essentially  ignore  it.  If  a  known  discrete 
falls  within  a  track  gate,  that  track  will  be  treated  as  if  in  a  shadow  and  will  not  be  updated. 

In  Figure  2-7,  the  same  target  speed,  update  time,  and  acceleration  noise  have  been  used 
as  in  the  shadow  rule  examples.  A  discrete  has  been  inserted  in  the  target  trajectory  as  shown. 
As  indicated,  the  tracker  preserves  the  state  value  of  4  as  the  predicted, gate  positions  passes 
through  the  discrete. 
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Figure  2-7:  Tracker  Response  for  Target  Passing  Near  a  Discrete 
Discrete  Rule  Applied 

The  previous  discussion  looked  at  the  implementation  of  four  specific  tracking  rules. 
However,  many  more  rules  that  were  considered  relevant  to  the  STAP  problem  were  developed 
during  the  course  of  this  tracking  study.  In  addition  to  the  rule  topics  discussed  here,  issues 
regarding  the  assignment  of  target  priority,  detection  threshold,  state  promotion  logic,  and  other 
features  were  considered.  A  collection  of  twenty-five  such  rules  are  presented  in  the 
“Knowledge-Based  Tracker  Rule  Book”  that  is  included  in  Volume  V  of  this  report.  Included 
with  each  rule  is  a  discussion  of  its  rationale,  its  impact  on  the  overall  knowledge  based  system 
and  interface  requirements  that  might  exist  between  the  tracker,  controller  and  radar. 

3.0  Software  Description 

All  tracking  software  used  in  this  effort  has  been  written  in  Matlab  5.1.  A  total  of  49 
modules  are  listed  and  documented  below.  These  can  be  subdivided  into  three  groupings 
consisting  of  user  interface,  main  tracker,  and  tracker  utilities.  The  first  two  modules  belong  to 
the  user  interface  group.  These  are  Tracker_GUI_7,  which  provides  the  graphical  user  interface 
for  the  overall  tracking  simulation,  and  Setup_Tracker_6,  which  sets  up  the  interface  between  the 
tracker  and  the  scenario  generator.  The  main  tracker  group  contains  the  multiple  target  tracking 
module,  Run_Tracker_8,  which  carries  out  all  prediction,  smoothing,  and  association  operations 
on  report  data  on  a  per  scan  basis.  Finally,  the  tracker  utilities  group  contains  all  of  the  support 
modules  that  are  used  by  the  tracker.  It  should  be  noted  that  some  of  the  modules  supporting  the 
user  interface  group  have  not  been  documented  since  that  portion  of  the  software  is  going  to  be 
removed  by  the  customer  and  replaced  with  their  own  hooks  into  the  STAP  simulation  software. 
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User  Interface  Modules 


Tracker_GUI_7 

Setup_Tracker_6 

Main  Tracker  Module 
Run_Tracker_8 

Tracker  Utility  Modules 

put_Table_R_2 

get_Table_R_2 

put_Table_Trk_4 

get_Table_Trk_4 

add_Table_Trk_4 

Cleanup_Table_Trk 

print_Table_Trk 

Display  _Trk_Data_5 

State_Update 

set_Trackstate 

set_TrackHM_stat 

set_Track_gates 

Shadow_Updates 

set_Track_Shadow 

predict_shlx_x 

Kal_b_pred 

Kal_c_pred 

smooth_shlx 

Kal_b_smth 

Kal_c_smth 

first_smooth_shlx 

abtrack_init 

track  init  2 
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Module  Name: 
Calling  Module: 
Called  Modules: 


Inputs: 

Outputs: 

Globals: 


Tracker_GUI_7 

None 

All_Defaults_4 

Demo_Defaults_4 

Read_Scenario_2 

W  rite_Scenario_2 

Track_Data_Gen_4 

Radar_menu_3 

Tracker_Opt_menu_2 

Setup_Tracker_6 

Run_Tracker_8 

None 

None 


*  *  *Radar_menu_3  *  *  * 
first_in_Radar 

f_Ghz  N_hits  SNR0_dB  Prf_kHz 

DRng_m  Az_mrad_3dB  Dop_hz 

Pfa  SNR_fa_dB 

***Scenario_menu*** 

script_name 

Tracker_Opt_resp 

***Track_Filter_menu*** 

Tfilt_resp 

*  *  *Tracker_param_menu*  *  * 
first_in_Tr_parm 
Tfilt_resp 

Mult_LP  Mult_HP  radius_TENT_LP  radius_TENT_HP 

Man_amt  alpha  beta  Pri_T 

Turn_int 

gate_case 

*  *  *Tracker_prom_menu  *  *  * 
first_in_Tr_prom 

Hit_Tbl_LP  Miss_Tbl_LP  Hit_Tbl_HP  Miss_TbLHP 
fig_no_T 

Time  X-meas  Y_meas  pDop_meas  report_ID 
X_tru  Y_trum  pDop_tm 
Sll  S12S13S22  S23  S33 
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ang_pred 

ang_trk 

gate_on 

print_Trk_Tbl 

text_on 

N_shadow 

x_sh_LO 

x_sh_HI 

y_sh_LO 

y_sh_HI 

In_Shadow 

Nshadow_dwens 


Angle  (radians)  of  prediction  point 
Angle  (radians)  of  track 
Draw  track  gate  flag  (1=>  yes) 

Print  track  gate  flag  (1=>  yes) 

Label  track  with  text  (1=>  yes) 

Number  of  shadow  zones  in  scenario 
Array  of  low  x  shadow  tile  values 
Array  of  high  x  shadow  tile  values 
Array  of  low  y  shadow  tile  values 
Array  of  high  y  shadow  tile  values 
In  shadow  flag  (1=>  yes) 

Number  of  consecutive  dwells  in  shadow 


Globals: 

Mult 

Table_R 

Table_T 

DROPPED  TENT  FIRM 
LPHP 

SUMl_ERROR  SUM2_ERROR 

P_DET 

HP_Buf 

dBuf 

Display_cnt 


Description: 

Initializes  and  sets  up  variables  to  be  used  by  the  Run_Tracker_8  program.  Puts  scenario 
generator  outputs  in  a  form  usable  by  tracker.  Also  adds  false  alarms  to  scenario 
generation  data. 
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Module  Name: 


Run  Tracker  8 


Calling  Modules:  None  (script  file) 

Called  Modules: 

put_Table_R_2 

get_Table_R_2 

get_Table_Trk_4 

put_Table_Trk_4 

add_Table_Trk_4 

predict_shlx_x 

first_smooth_shlx 

smooth_shlx 

Get_Semiaxes_3 

Test_E_Gate 

Test_Smile_Gate 

get_def_Cov_manu 

Assign_Cov_manu 

State_Update 

set_Trackstate 

set_TrackHM_stat 

Pred_Shadow_Test 

Shadow_Update 

Discrete_Test 

set_Track_Shadow 

get_track_ang 

get_pred_ang 

update_error 

TR_Assoc_Max_T 

Display_Trk_Data_5 

print_Table_T 

Load_dBuf_4 

plot_dBuf 

Cleanup_Table_Trk 

Inputs:  setup  by  Tracker_GUI_7  and 

Setup_Tracker_6 

Outputs:  None 

Globals:  None 

Description: 

The  multiple  target  tracker  processes  report  data  collected  during  each  scan  and  performs 
three  basic  functions:  track-report  correlation,  association,  and  maintenance. 
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Module  Name: 


Calling  Module: 

Called  Modules: 

Inputs: 

Time 

X_meas,  Y_meas 
FDop_meas 
X_tru,  Y_tru 
FDop_tru 

Sll . S33 

sig_Rng_km 

sig_Az_rad 

sig_Rdot_kmps 

freq_GHz 

Det_Level 

Priorjntel 

Manu_intel 

report_ID 

Nrep 

Called  Modules: 

Inputs: 


Time 

X_meas,  Y_meas 

FDop_meas 

X_tru,  Y_tru 

FDop_tru 

Sll,  ...,S33 

sig_Rng_km 

sig_Az_rad 

sig_Rdot_kmps 

freq_GHz 

Det_Level 

Prior_intel 

Manu_intel 

report_ID 

Nrep 

scan 

Outputs: 


put_Table_R_2 

Run_Tracker 

None 


Time  of  radar  blips 
Measured  position  (km) 

Measured  pseudo  Doppler  (km*km/sec) 

True  position  (km) 

Truer  pseudoDoppler  (km*km/sec) 

Measurement  Covariance 
Range  Measurement  error  (km) 

Azimuth  measurement  error  (radians) 

Range  rate  measurement  error  (km/sec) 

Radio  frequency  of  radar  (gHz) 

Detection  level 
Priority  status 
Maneuver  status 
Trajectory  ID  of  report 

Number  of  reports  calling  Module:  Run_Tracker 
None 


Time  of  radar  blips  (sec) 

Measured  position  (km) 

Measured  pseudo  Doppler  (km*km/sec) 

True  position  (km) 

True  pseudoDoppler  (km*km/sec) 

Measurement  Covariance 
Range  Measurement  error  (km) 

Azimuth  measurement  error  (radians) 

Range  rate  measurement  error  (km/sec) 

Radio  frequency  of  radar  (gHz) 

Detection  level 
Priority  status 
Maneuver  status 
Trajectory  ED  of  report 

Number  of  reports  calling  Module:  Run_Tracker 
Current  scan  index 

Table_R 
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Globals: 


Table  R 


Description: 

Load  the  report  table  buffer  for  the  current  scan.  All  arrays  are  indexed  as  (report,  scan) 
and  were  created  by  a  scenario  generator. 
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Module  Name: 


get_TabIe_R_2 
Calling  Module:  Run_Tracker_8 

Called  Modules:  None 

Inputs: 

rep  report  index 

Outputs: 

T_m 

X_m,  Y_m 
FDop_m 
Cov_m_vec 
sig_rng_km_m 
sig_cmg_km_m 
sig_Rdot_kmps_m 
freq_GHz 
Det_m 
priority_in 
maneuver_in 
report_source 
X_tru_m,  Y_tru_m 
FDop_tru_m 
error_status 

Globals:  Table_R 

Description: 

Fetch  report  table  data  for  current  scan. 


Blip  time  (sec) 

Measured  position  (km) 

Measured  Doppler  (km*km/sec) 
Measured  covariance  vector 
Range  measurement  error  (km) 

Cross  range  measurement  error  (km) 
Range  rate  measurement  error  (km?sec) 
Radar  rf  frequency  (gHz) 

Detection  level 
Priority  Status 
Maneuver  status 
Trajectory  ID 
True  position  (km) 

True  pseudoDoppler  (km*km/sec) 
not  used 
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Module  Name: 

put_Table_Trk_4 

Calling  Module: 

Run_Tracker_8 

Called  Modules: 

None 

Inputs: 

trk 

Track  index 

TID 

Track  ID  number 

status 

Track  status  (drop,  tentative,  firm) 

state 

Track  state  (quality  index;  0  =>  dropped) 

report 

Captured  report  index 

scanl 

Scan  index  when  track  became  firm 

trk_type 

Tracker  type  (alpha  beta,  unc/extd  Kalman) 

Z_m 

Measurement  state  vector 

Cov_m_vec 

Measurement  covariance  vector 

T_m 

Measurement  of  time  of  captured  report 

priority 

Priority  status  of  track 

maneuver 

Maneuver  status  of  track 

Z_p 

Prediction  state  vector 

Cov_p_vec 

Prediction  covariance  vector 

Z_s 

Smoothed  state  vector 

Cov_s_vec 

Smoothed  covariance  vector 

HM_flg 

Track  capture  flag  (0  =>  miss,  1  +.  hit) 

ang_trk 

Track  angle  (radians) 

x_s_2LST 

Last  track  x  position  (km) 

y_s_2LST 

Last  track  y  position  (km) 

T_s_2LST 

Last  track  time  (sec) 

semi_rng_T 

along  range  semi  axis  (km) 

semi_crng_T 

cross  range  semi  axis  (km) 

semi_trk_T 

along  track  semi  axis  (km) 

semi_ctrk_T 

cross  track  semi  axis  (km) 

In-Shadow 

Shadow  status  flag  (0  =>  not  in) 

Nshadow_dwells 

Number  of  successive  dwells  in  shadow 

gate_case 

Track  gating  choice  (0:3) 

Accel_max 

Maximum  acceleration  (km/sec*sec) 

sig_rng_km 

range  measurement  error  (km) 

sig_cmg_km 

cross  range  measurement  error  (km) 

sig_Rdot_kmps2 

range  rate  meas  error  (km*km/sec) 

Outputs: 

error_status 

Table_T 
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Globals: 


Table_T 


Description: 

Updates  track  table  data  for  current  scan  and  report  index. 
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Module  Name: 
Calling  Module: 
Called  Modules: 
Inputs: 

Outputs: 

TID 


get_Table_Trk_4 

Run_Tracker_8 

None 

trk 

same  as  inputs  to  put_Table_Trk_4 
same  as  inputs  to  put_Table_Trk_4 


sig_Rdot_kmps2 

Globals:  Table_T 

Description: 

Fetch  track  table  data  for  current  scan  and  report  index. 
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Module  Name: 


add_TabIe_Trk_4 
Calling  Module:  Run_Tracker_8 

Called  Modules:  None 

Inputs: 

Ntrk 

num_new_trks 

TID_New  same  as  inputs  to  put_Table_Trk4 

sig_Rdot_kmps2 

Outputs: 

error_status 

Table_T 

Globals  Table.T 

Description: 

Insert  new  tentative  tracks  into  track  table. 
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Module  Name: 
Calling  Module: 
Called  Modules: 
Inputs: 

Outputs: 

Ntrk 

Globals: 

Description: 


Cleanup_Table_Trk 

Run_Tracker_8 

None 

None 

Number  of  valid  tracks 
Table_T 


Updates  track  table  for  next  scan.  Sorts  Table_T  by  column  3  (state)  and  removes 
zero  states.  Counts  number  of  remaining  valid  tracks  and  returns  this  value. 


Module  Name: 


Calling  Module: 
Called  Modules: 
Inputs: 

Ntrk 

sup 

Outputs: 


print_Table_T 

Run_Tracker_8 

None 

Number  of  track  sets  to  be  printed  (max  =10) 
vector  of  Table_T  indices  to  be  printed 

prints  to  command  tool 


Globals:  Table_T 

Description: 

Prints  track  table  information  as  columns  (one  per  track).  Support  vector  sup  is  subset  of 
{1:86}  corresponding  to  fields  in  Table_T. 
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Module  Name: 
Calling  Module: 
Called  Modules: 


Inputs: 

Ntrk 

scan 

scan_period 

sig_rng_km_nom 

sig_Az_rad_nom 

plot_vec 

x_v,  y_v 

dx,  dy 

gate_on 

text_on 

fig_no 

symboLl 

symbol_2 

semi-max 

Outputs: 


Display_Trk_Data_5 

Run_Tracker_8 

Draw_Grid 

get_Table_Trk_4 

get_Table_R_2 

Assign_Cov_manu 

predict_shlx_x 

get_sig_track 

get_pred_ang 

E_Gate_Semiaxes 

Draw_Smile_Gate 

Draw_E_Gate_2 


Number  of  tracks  to  be  plotted 

current  scan  index 

scan  time  (sec) 

not  used 

not  used 

plot  scale:  (xmin,  xmax,  ymin,  ymax) 
X  and  y  grid  point  sets 
X  and  y  grid  spacing 
Draw  track  gate  flag  (1  =>  yes) 

Print  Text  gate 
plot  figure  number 
rng-cross  rng  aligned  gate  symbol 
trk-cross  trk  aligned  gate  symbol 
maximurri  gate  size  allowed 

graphical  display  in  figure  (fig_no) 


Globals:  Display_cnt 

Table_R  Table.T 

DROPPED  TENT  FIRM  LP  HP 

Mult_LP  Mult.HP  radius_TENT_LP  radius_TENT_HP 

sig_trk_km  sig_ctrk_km  gate_case 

Cov_man20  Cov_man30  Cov_man_acc_x  Cov_man_acc_y 


Description: 

The  module  serves  as  a  diagnostic  and  display  tool  for  the  multi  target  tracker.  It 
displays  predicted  and  measured  positions  prior  to  track  table  cleanup  and  draws  track 
gates  centered  on  predicted  positions. 
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Module  Name:  set_Trackstate 

Calling  Module:  Run_Tracker_8 

Called  Modules:  None 

Inputs: 

state 

trk 

Outputs:  None 

Globals:  Table_T 

Description: 

Inputs  state  value  into  track  table  at  index  trk. 
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Module  Name: 


set_TrackHM_stat 
Calling  Module:  Run_Tracker_8 

Called  Modules:  None 

Inputs: 

HM_flg  Capture/shadow  status  of  track 

trk  Track  index 

Outputs:  None 

Globals:  Table_T 

Description: 

Inputs  the  track  capture/shadow  status  into  track  table. 
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Module  Name: 


set_Track_gates 
Calling  Module:  Run_Tracker_8 

Called  Modules:  None 

Inputs: 

semi_rag_T  Along  range  semi  axis  of  ellipse  gate 

semi_crng_T  Cross  range  semi  axis 

semi_trk_T  Along  track  semi  axis 

semi_ctrk_T  Cross  track  semi  axis 

trk  Track  table  index 

Outputs:  None 

Globals:  Table.T 

Description: 

Inputs  ellipse  gate  semi  axes  into  track  table  corresponding  to  index  trk. 


Module  Name: 


Shadow_Update 
Calling  Module:  Run_Tracker_8 

Called  Modules: 

set_Track_Shadow 

set_TrackHM_stat 

Inputs: 

trk  Track  index  into  Table_T 

In_Shadow  Shadow  status  flag  (1  =>  in  shadow) 

Nshadow_dwells_last  Number  of  scans  input  trk  in  shadow 
HM_flg_last  Capture/shadow  status  of  input  trk 

Outputs:  None 

Globals:  Table.T 

Description: 

Updates  track  shadow  status  corresponding  to  index  trk.  If  track  is  in  a  shadow,  the 
number  of  shadow  dwells  is  incremented  and  the  capture  status  is  set  to  the  number  2. 
Otherwise,  the  number  of  dwells  is  et  to  zero  and  the  capture  status  maintained. 
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Module  Name: 


set_Track_Shadow 
Calling  Module:  Run_Tracker_8  Shadow_Update 

Called  Modules:  None 

Inputs: 


Nshadow_dwells  Number  of  consecutive  dwells  in  shadow 
In_Shadow  in  shadow  status  flag 

trk  Index  into  Table_T 

Outputs:  None 

Globals:  Table_T 

Description: 

Inputs  shadow  parameters  into  track  table  corresponding  to  index  trk.  These  include  (1) 
the  #  of  consecutive  dwells  of  track  in  a  shadow,  and  (2)  the  in-shadow  status  flag  (0  => 
not  in  shadow,  1  =>  in  shadow). 
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Module  Name: 
Calling  Module: 
Called  Modules: 


Inputs: 

type 

Z_s_last 
Cov_s_last_vec 
Cov_man_x 
Cov_man_y 
Cov_man 
dT 

Outputs: 

Z-P 

Cov_p_vec 
Globals: 

Description: 

Computes  predicted  state  vector  and  Covariance  array,  Z_p  and  Cov_p_vec,  respectively. 
The  method  used  depends  on  the  tracking  filter  type  specified.  The  maneuver  covariance 
matrices  provide  a  means  for  inputting  acceleration  noise  to  increase  track  gate  size  for 
the  Kalman  filter  cases 

Equations: 

alpha  beta: 

X_p  =  x_s_last  +  vx_s_last  *  dT 
Y_p  =  y_s_last  +  vy_s_last  *  dT 
Vx_p  =  vx_s_last 
Vy_p  =  vy_s_last 

where  the  components  X_p,  Vx_p,  etc.,  are  related  to  the  state  vectors  Z_p,  Z_s_last,  etc. 
by  stacking  the  x  and  y  components  of  position,  velocity,  and  acceleration 

Z=[X;Vx;Ax;Y;Vy;Ay] 


predict_shlx_x 

Run_Tracker_8 


Kal_b_pred 

Kal_c_pred 


Tracking  filter  type  (1:3  +.  alpha  beta,  uncoupled  Kalman,  and 
extended  Kalman) 

Smoothed  state  vector  at  last  scan 
Smoothed  covariance  array  at  last  scan 
Maneuver  covariance  matrix  of  x  component 
Maneuver  covariance  matrix  of  y  component 
Maneuver  covariance  matrix  of  extended  Kal 
Time  interval  from  last  update  to  present 


Prediction  state  vector 
Prediction  covariance  matrix 

None 
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Uncoupled  Kalman: 

see  Kal_b_pred 

Extended  Kalman: 

see  Kal_c_pred 
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Module  Name: 


Kal_b_pred 

Calling  Module:  predict_shlx_x 

Called  Modules:  None 

Inputs: 

Z_s_in 
dt_in 

Cov_manu 
Cov_s_in 


Outputs: 

Z_p_out 
Cov_p 

Globals: 

Description: 

Computes  3x1  prediction  state  vector  and  3x3  prediction  covariance  matrix  for  the  case 
of  an  uncoupled  Kalman  filter.  All  input  and  output  state  vectors  are  assumed  to  be  3  x  1, 
and  all  covariance  matrices  are  3  x  3. 

Equations: 

1  dT  0.5{dTf' 

0  1  dT  ;  3  X  3  state  transition  matrix 

0  0  1 

Z_p_out=<t)*Z_s_in 
Cov_p=(|)*Cov_s_in*(l)’+Cov_manu 
((])’  denotes  the  transpose  of 


3x1  Smoothed  input  state  vector  (either  component) 
[Position;  Velocity;  Acceleration] 

Time  since  last  update 

3x3  Covariance  matrix  of  maneuver  noise 
3x3  Smoothed  input  covariance  matrix  (either  component) 
[Position;  Velocity;  Acceleration] 


3x1  Prediction  state  vector  (corresponding  component) 
3x3  Prediction  covariance  matrix 

Table_T 
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Module  Name: 


KaLc_pred 

Calling  Modules;  predict_shlx_x 

Called  Modules:  None 

Inputs: 


Z_s_in 

4x1 

Smoothed  input  state  vector 
[X_s_in;  Vx_s_in;  Vy_s_in] 

dt_in 

Time  since  last  update 

Cov_manu 

4x4 

Covariance  matrix  of  maneuver  noise 

Cov_S_in 

4x4 

Smoothed  input  covariance  matrix 

Outputs: 


Z_p_out 

4x1 

Prediction  state  vector 

[X_p_out;  Vx_p_out;  Y_p_out;  Vy_p_out 

Cov_p 

4x4 

Prediction  covariance  matrix 

Globals;  None 

Description: 

Computes  4x1  prediction  state  vector  and  4x4  prediction  covariance  matrix  for  the  case 
of  an  extended  Kalman  filter.  All  input  and  output  state  vectors  are  assumed  to  be  4  x  1, 
and  all  covariance  matrices  are  4  x  4. 


(([)’ denotes  the  transpose  of  (j)) 
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Module  Name: 
Calling  Module: 
Called  Modules: 


Inputs: 

type 

Z_m 

Z_p 

Cov_p_vec 

Cov_m_vec 

dT 

Outputs: 

Z_s 

Cov_s_vec 

Gain_mat 

Globals: 

Description: 


smooth_shlx 

Run_Tracker_8 


Kal_b_smth 

Kal_c_smth 


Tracking  filter  type  (1:3  =>  alpha  beta,  uncoupled  Kalman,  and 
extended  Kalman) 

Measurement  vector  (x;y;pDop) 

Prediction  state  vector  (X;Vx;Ax;Y;Vy;Ay) 

Prediction  cov  array  (Cov_p_x;Cov_p_y)reshaped  to  1  x  18 

Meas  cov  array  1x9 

Time  interval  from  last  update  to  present 


Smoothed  state  vector  (X;Vx;Ax;Y;Vy;Ay) 

Smoothed  cov  array  (Cov_s_x;Cov_s_y)reshaped  to  1  x  18 
Gain  matrix  (not  used  by  tracker) 

alpha  beta 


Computes  smoothed  state  vector  and  Covariance  array,  Z_s  and  Cov_s_vec,  respectively. 
Also  returns  a  gain  matrix  (not  normally  used)  which  can  be  used  for  diagnostic 
purposes.  The  smoothing  method  used  depends  on  the  tracking  filter  type  specified. 


Equations: 

alpha  beta: 


where 


Dx  =  x_m  -  X_p 

x_s  =  X_p  +  alpha  *  Dx 

Vx_s  =  Vx_p  +  beta  *  Dx/dT 


Dy  =  y_m  -  Y_p 

y_s  =  Y_p  +  alpha  *  Dy 

Vy_s  =  Vy_p  +  beta  *  Dy/dT 

Z_m  =  [x_m;  y_m;  pDop_m) 


and  Z_p  and  Z_s  are  6  x  1  vectors  of  the  position,  velocity  and  acceleration  components 
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Uncoupled  Kalman:  see  Kal_b_smth 


Extended  Kalman:  see  Kal_c_smth 
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Module  Name: 

Kal_b_smth 

Calling  Module: 

smooth_shlx 

Called  Modules: 

None 

Inputs: 

Z_m 

3x1 

Measurement  vector  [x_m  ;y_m  ;pDop_m] 

Z_p_out 

3x1 

Prediction  state  vector  (either  component) 
[Position;  Velocity;  Acceleration] 

Cov_p 

3x3 

Covariance  matrix  of  maneuver  noise 

Cov_s_in 

3x3 

Prediction  covariance  matrix 

var_meas 

Measurement  variance 

Outputs: 

Z_s_out 

3  X  1 

Smoothed  output  state  vector  (either  component) 
[Position;  Velocity;  Acceleration] 

Cov_s_out 

3x3 

Smoothed  output  covariance  matrix 

Res 

3x1 

Residual  error  vector  (Innovations) 

dist  Res  2 

Statistical  distance 

S 

3x1 

Gain  matrix 

Globals: 

None 

Description: 

Computes  3  x  1  smoothed  state  vector  and  3  x  3  smoothed  covariance  matrix  for  the  case 
of  an  uncoupled  Kalman  filter.  All  input  and  output  prediction  and  smoothed  state 
vectors  are  assumed  to  be  3  x  1,  and  all  covariance  matrices  are  3  x  3. 

Equations: 

(for  either  x  or  y  component) 

M  =  [1  0  0];  Measurement  Matrix  (1  Measurement  x  3  states) 

I_Cov_Res  =  1/(M  *  Cov_p  *  M’  +  var_meas) 

S  =  Cov_p  *  M’  *  I_Cov_Res;  Gain  Matrix  (3  states  x  1  measurement) 

Z_s_out  =  Z_p_out  +  S  *  (Z_m  -  M  *  Z_p_out) 

Cov_s_out  =  (I  -  S  *  M)  *  Cov_p 

Res  =  (Z_m  -  M  *  Z_p) ;  Innovations 

dist_Res_2  =  Res  *  I_Cov_Res  *  Res';  Statistical  distance 
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Module  Name: 


Kal_c_smth 
Calling  Modules;  smooth.shlx 

Called  Modules:  None 

Inputs: 


Z_m 

3x1 

Measurement  vector  (3x1) 
[x_m;y_m;pDop_m] 

Z_p_out 

4x1 

Prediction  state  vector  [X;  Vx;  Y;  Vy] 

Cov_p 

4x4 

Prediction  covariance  matrix 

Cov_meas 

3x3 

Measurement  Covariance 

Outputs: 


Z_s_out 

4x1 

Smoothed  output  state  vector  [X;  Vx;  Y;  VyJ 

Cov_s_out 

4x4 

Smoothed  output  covariance  matrix 

Res 

3x1 

Residual  error  vector  (Innovations) 

dist_Res_2 

Statistical  distance 

S 

4x3 

Gain  Matrix 

Globals:  None 

Description: 

Computes  4  x  1  smoothed  state  vector  and  4  x  4  smoothed  covariance  matrk  for  the  case 
of  an  extended  Kalman  filter.  All  input  and  output  prediction  and  smoothed  state  vectors 
are  assumed  to  be  4  x  1,  and  all  covariance  matrices  are  4  x  4. 


Equations: 


M  = 


1 

0 


0  0  0 

0  1  0 

^32  ^33  ^34 


Measurement  matrix  (3  Measurements  x  4  states) 


where 

M3i=0.5*Z_p_out  (2) 
M32  =  0.5  *  Z_p_out  (1) 
M33  =  0.5  *  Z_p_out  (4) 
M34  =  0.5  *  Z_p_out  (3) 
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Note  that 

(pDop)est  =  M31  *  X_p  +  M32  *  Vx_p  +  M33  *  Y_p  +  M34  *  Vy_p 
I_Cov_Res  =  (M  *  Cov_p  *  M’+  Cov_meas)-l 

S  =  Cov_p  *  M’  *  I_Cov_Res ;  4  x  3  Gain  Matrix  (4  states  x  3  measurements) 

Res  =  (Z_m  -  M  *  Z_p_out) ;  3  x  1  Innovations 

Z_s_out  =  Z_p_out  +  S  *  Res  4  x  1 

Cov_s_out  =  (I  -  S  *  M)  *  Cov_p  4x4 

dist_Res_2  =  Res;  *  (I_Cov_Res)  *  Res 
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Module  Name: 


first  smooth  shlx 


Calling  Module:  Run_Tracker_8 

Called  Modules: 

abtrack_init 

track_init_2 

track_init_c 

Inputs: 

type  Track  filter  type  (1:3  =>  alpha  beta,  Uncoupled  Kalman,  and 

extended  Kalman 

Z_m_last  Measured  vector  (last  update) 

Z_m  Measured  vector  (current) 

Cov_m_vec  Meas  Covariance  array  (current) 

dT  Time  interval  (sec)  from  last  update 

Outputs: 

Z_s_out  Smoothed  state  vector 

Cov_s_vec  Smoothed  covariance  array 

Globals:  None 

Description: 

Performs  smoothing  of  tentative  tracks 
Equations: 

alpha  beta: 

see  abtrack_init 
Uncoupled  Kalman:  see  track_init_2 

Extended  Kalman:  see  track_init_c 
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Module  Name: 


abtrack  init 


Calling  Module:  first_smooth_shlx 

Called  Modules :  N one 

Inputs: 

(xl,yl) 

(x2,  y2) 
dT12 

Outputs: 

(x_s,  y_s)  Smoothed  track  position 

(vx_s,  vy_s)  Smoothed  tracked  velocity 

ang_trk  Track  angle  (radians) 

Globals:  None 

Description: 

Does  a  two  point  initialization  of  an  alpha  beta  tracker  using  current  and  last  track 
positions  and  their  time  interval. 

Equations: 

Dx  =  (x2-xl) 
x_s  =  x2 
vx_s  =  Dx/dT12 

Dy  =  (y2  -  yl) 
y_s  =  y2 
vy_s  =  Dy/dT12 

ang_trk  =  arctan  (Dy/Dx) 


Last  updated  track  position 
Current  track  position 
Time  interval  (sec) 
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Module  Name: 
Calling  Module: 
Called  Modules: 
Inputs: 


Delt 

z_pos_l 

z_pos_2 

var_pos_m_2 

Outputs: 

Z_s_in 

P_s_in 

Globals: 


track_init_2 

first_smooth_shlx 

None 


Time  (sec)  between  current  and  last  update 
Last  position  component  update  (x  or  y) 
Current  position  component  (x  or  y) 
Measurement  variance  of  position  component 


Smoothed  state  vector  (3x1) 
Smoothed  covariance  matrix  (3  x  3) 

None 


Description: 

Does  a  two  point  initialization  of  an  uncoupled  Kalman  tracker.  This  routine  is  applied 
separately  to  both  the  x  and  y  components  of  target  motion. 


Equations: 


pos_s  =  z_pos_2 

vel_s  =  (z_pos_2  -  z_pos_l)/Delt 
Z_s_in  =  [pos_s;  vel_s;  0] 


P 


s  in  = 


PXy 

Px 

0 


12 


Px„  0 
PX22  0 
0  0 


where 

Pxii  =  var_pos_m_2 
Pxi2  =  Pxii/Delt 
Px22-Pxii/(Delt)2 
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Module  Name: 


track_init_c 


Calling  Module: 
Called  Modules: 
Inputs: 


Delt 
(xl,yl) 
(x2,  Y2) 
sig_xx 

sig-yy 

sig_xy 

Outputs: 

Z_s_in 

P_s_in 


first_smooth_shlx 

None 


Time  (sec)  between  current  and  last  update 
Last  position  update 
Current  position 

Measurement  covariance  xx  comp 
Measurement  covariance  yy  comp 
Measurement  covariance  xy  comp 


4x1  Smoothed  state  vector 
4x4  Smoothed  covariance  matrix 


Globals:  None 

Description: 


Does  a  two  point  initialization  to  be  used  with  the  extended  Kalman  track  filter. 


Equations: 

Let  T  =  Delt,  T2  =  T^  _  _ 

sig  _x  =  ^sig_  xx,sig  _y  =  ylsig_yy 


pos_sx  =  x2,  pos_sy  =  y2;  smoothed  positions 

veLsx  =  (x2  -  xl)/T,  vel_sy  =  (y2  -  yl)/T  ;  smoothed  velocities 

Z_s_in  =  [pos_sx;  vel_sx;  pos_sy;  vel_sy] 


^12 

^13 

Pl4 

in  = 

^13 

P22 

P2J 

^24 

^23 

^33 

^4 

Pl4 

P24 

P34 

P44 

where 

Pii  =  sig_xx,  Pi2  =  1.5  *  Pii/T,  Pi3  =  sig_xy;  Pu  =  1.5  *  sig_x  *  sig_y/T 
P22  =  6.5  Pii/T2;  P23  =  Pi4.  P24  =  sig_xy/T2 
P33  =  sig_yy,  P34  =  1.5  *  P33/T 
P44  =  6.5  *  P33A’2 
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Module  Name: 
Calling  Module: 
Called  Modules: 
Inputs: 


trk_type 

scan 

manu_rep 

rep_source 

Outputs: 

Cov_man 

Cov_man_x 

Cov_man_y 

Globals: 


Assign_Cov_inanu 

Run_Tracker_8 

get_def_Cov_manu 


Track  filter  type 
Current  scan 

Maneuver  anticipation  status  flag 
Maneuver  anticipation  interval  (scan  lo,  scan  hi) 


Maneuver  covariance  matrix 


Cov_man20 

Cov_man30 

Cov_man_acc_x 

Cov_man_acc_y 

Turn_int 


Description: 

Sets  up  the  maneuver  covariance  matrix.  If  a  maneuver  has  been  anticipated  for  the 
current  scan  then  maneuver  noise  is  applied.  Otherwise  a  small  value  (essentially  zero)  is 
loaded  into  the  matrix. 


Equations: 

Accel_max  =  Maximum  maneuver  acceleration  (km/s ) 

var_acc  =  (Accel_max)^ 

var_vel  =  (Accel_max  *  scan_period)^ 

var_pos  =  (0.5  *  Accel_max  *  (scan_period)^)^ 


Cov_man  = 


y2Lr_pos 

0 

0 


0 

var_  vel 
0 


0 

0 

\BT_acc 


trk_type  ^  3 


46 


Cov_man_x  =  Cov_man_y  =  Cov_man 
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Module  Name: 

Calling  Module: 

Called  Modules: 

Inputs: 

trk_type  Track  filter  type 

Outputs: 

Cov_man  Maneuver  covariance  matrix 

Cov_man_x 

Cov_man_y 

Globals: 

Cov_man20 

Cov_man30 

Description: 

Loads  default  values  into  maneuver  covariance  matrices. 
Equations: 

Cov_man  =  Cov_man20;  trk_type  3 
Cov_man  =  Cov_man30;  trk_type  =  3 
Cov_man_x  =  Cov_man_y  =  Cov_man; 


get_def_Cov_manu 

Assign_Cov_manu 

None 


where 


Cov_man20  =  le-10  *  eye  (3) 
Cov_man30  =  le-6  *  eye  (4) 
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Module  Name: 


Get  Semiaxes  3 


Calling  Module: 


Called  Modules: 


Inputs: 


In_Shadow 

In_Shadow_last 

N  shadow_dwells_ 

gate_case 

trk_type 

Mult 

Accel_max 

dT 

ang_pred 

ang_trk 

sig_rng_km_m 

sig_crag_km_m 

Z_s_last 

Z_p 

Cov_p_vec 
senii-rng_T_last 
semi_crng_T_las  t 
semi_trk_T_last 
semi  ctrk  T  last 


semi  max 


Outputs: 


semi_rng_T 
semi_crng_T 
semi_trk_T 
semi  ctrk  T 


Globals: 


Description: 


Run  Tracker  8 


E_Gate_S  emiaxes_3 

get_sig_track 

Scale_Shadow_gates 


Current  in  shadow  flag  (1  =>  yes) 

Last  update  in  shadow  flag 

last  Number  of  successive  dwells  in  shadow  as  of  last  update 
Gate  type  flag  (0:3) 

Track  filter  type  ( 1 :3) 

Gate  size  multiplier 

Maximum  assumed  acceleration 

Time  (sec)  since  last  update 

Angle  (rad)  of  predicted  point  wrt  radar  cs 

Angle  (rad)  of  track  wrt  radar  cs 

Range  measurement  error  (km) 

Cross  range  measurement  error  (km) 

Smoothed  state  vector  (last  update) 

Prediction  state  vector  (current) 
prediction  covariance  array  (current) 

Along  range  semi  axes  (last) 

Cross  range  semi  axes  (last) 

Along  track  semi  axes  (last) 

Cross  track  semi  axes  (last) 

Maximum  allowed  semi  axis  size 


Along  range  semi  axes  (current) 
Cross  range  semi  axes  (current) 
Along  track  semi  axes  (current) 
Cross  track  semi  axes  (current) 

None 


Computes  the  semi  axes  for  the  two  orientations  of  elliptical  gates  used.  The  first  type  is 
oriented  along  and  across  range  and  the  second  is  oriented  along  and  across  track.  Tracks 
that  are  in  a  shadow  have  their  gate  sizes  frozen.  Otherwise,  track  gates  sizes  are 
determined  as  a  function  of  the  gate  case  and  track  filter  type  selected.  TTiis  function  is 
done  by  the  routine  E_Gate_Semiaxes_3. 
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Module  Name: 

E_Gate_Seimaxes_3 

Calling  Module: 

Get_Semiaxes_3 

Called  Modules: 

Rotate_xy2xpyp 

Inputs: 

gate_case 

trk_type 

Mult 

Z_p 

Cov_p_vec 

ang_pred 

sig_rng 

sig_cmg 

ang_trk 

sig_trk 

sig_ctrk 

semi_max 

Gate  type  flag  (0:3) 

Track  filter  type  (1:3) 

Gate  size  multiplier 

Prediction  state  vector  (current) 
prediction  covariance  array  (current) 
Angle  (rad)  of  predicted  point  wrt  radar  cs 
Range  measurement  error  (km) 

Cross  range  measurement  error  (km) 
Angle  (rad)  of  track  wrt  radar  cs 

Along  track  measurement  error  (km) 

Cross  track  measurement  error  (km) 
Maximum  allowed  semi  axes  size 

Outputs: 

semi_rng 

semi_crng 

semi_trk 

semi_ctrk 

Along  range  semi  axes  (km) 

Cross  range  semi  axes  (km) 

Along  track  semi  axes  (km) 

Cross  track  semi  axes  (km) 

Globals; 

None 

Description; 

Computes  semi  axes  for  elliptical  gates  oriented  along/cross  range,  and  oriented 
along/cross  track.  Results  depend  on  which  of  three  gate  cases  are  chosen  (cases  1  and  2 
require  a  Kalman  filter).  In  each  case  the  along/cross  track  gate  sizes  are  computed  as 
scaled  versions  of  the  along/cross  track  measurement  errors.  However,  the  along/cross 
range  oriented  gate  sizes  are  case  dependent.  Gate  case  1  combines  in  an  rss  fashion  (1) 
range  and  cross  range  measurement  errors  from  the  last  track  update,  and,  (2)  current 
prediction  covariance  estimates  in  x  and  y,  projected  onto  the  range/cross  range  axes. 
The  composite  gate  is  formed  by  multiplying  this  result  by  a  scale  factor.  Gate  case  2 
uses  only  the  scaled  prediction  covariance  estimates  in  x  and  y  to  forai  the  gate.  Gate 
case  3  uses  the  scaled  along/cross  range  measurement  errors  to  form  this  gate. 


Equations: 
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The  quantities  sig_mg  and  sig_crng  are  measurement  errors  along  range  and  cross  range 
computed  from  the  last  report  captured  by  a  given  track  and  stored  in  the  track  table.  The 
quantities  sig_trk  and  sig_ctrk  are  based  on  kinematical  assumptions  and  are  computed  in 
get_sig_track. 


alpha  beta  filter: 


semi_mg 

sig_mg 

semLcmg 

sig_cmg 

semi_trk 

sig_trk 

semi_ctrk 

_sig_ctrk_ 

1.5MMfr ,  (gate_case  <  3) 
Mult  ,  (gate_case  =  3) 


Kalman  filters: 


Get  prediction  uncertainty  along  x  and  y  directions  (Opx,  (Jpy)  from  the  covariance  array 
Cov_p_vec. 


(gate_case  =1) 

A0  =  (ang_pred  -  ang_trk) 

Compute  projections  of  Opx  and  Opy  onto  the  range-cross  range  axes. 


cos(A0) 

sinfA^) 

^px 

’  -sin(A^) 

cos(A^)_ 

Form  the  “root  sum  square”  (Rss)  of  Opx’  and  Opy-with  the  measurement  errors 

semi_mg  =  Mult  *  yj{apx'f  -i-(sig_mg)^ 

semLcrng  =  Mult  *  ^apx'f  +  (sig_cmg)^ 

semLtrk  =  Mult  *  sig_trk 
semi_ctrk  =  Mult  *  sig_ctrk 
(gate_case  =  2) 

6  =  ang_trk 

Compute  projections  of  Opx  and  OpyOnto  track-cross  track  axes. 


_r  cos(^) 

sinf^) 

~  -sin(0) 

cos(^)_ 

P  py . 
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semi_rng 

semi_cmg 

semi_trk 

=  Mult  * 

semi_ctrk 

(gate_case  =  3) 


semi_mg 

semi_cmg 

semi_trk 

semi_ctrk 


=  Mult  * 


sig_rng 

sig_cmg 


sig_rng 

sig_crng 

sig_trk 

sig_ctrk 
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Module  Name: 

Calling  Module: 
Called  Modules: 
Inputs: 


Z_s 

Accel_max 

dT 

Outputs: 

sig_trk_kni 

sig_ctrk_km 

Globals: 


get_sig_track 

Get_Semiaxes_3 

None 


Smoothed  state  vector 

Maximum  allowed  acceleration  (km/s*s) 

Time  (sec)  since  last  update 


Along  track  uncertainty  (km) 
Cross  track  uncertainty  (km) 

None 


Description: 

Computes  along/cross  track  uncertainties.  The  along  track  uncertainty  is  computed  as 
0.5*A*(dT)^  with  the  along  track  acceleration  A  assumed  bounded  by  0.5g.  The  across 
track  uncertainty  is  assumed  to  be  due  to  constant  speed  turning  only. 

Equations: 

Define  the  following  quantities: 

Acc_at  =  along  track  acceleration  <  Accel_max 
V_kmps  =  estimated  track  speed 
pkm  =  radius  of  curvature  of  turn  (km) 

Ang  =  angle  of  turn  voer  time  dT  (radians) 

pkm  =  (V_kmps)  /Accel_max 
ang  =  V_kmps  *  dT/pkm 

sig_trk_km  =  0.5  Acc_at  (dT)^ 
sig_ctrk_km  =  pkm  (1  -  cos  (Ang)) 
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Module  Name: 
Calling  Module: 
Called  Modules: 
Inputs: 


gate_scale 

semi-rng_T_in 

semi_cmg_T_in 

semi_trk_T_in 

semi_ctrk_T_in 

Outputs: 

semi_mg_T_out 

semi_crng_T_out 

semi_trk_T_out 

semi_ctrk_T_out 

Globals: 

Description: 

Multiplies  input  gate 


Scale_Shadow_Gates 

Get_Semiaxes_3 

None 


Gate  size  multiplier 
Along  range  semi  axis  (input) 
Cross  range  semi  axis  (input) 
Along  track  semi  axis  (input) 
Cross  track  semi  axis  (input) 


Along  range  semi  axis  (output) 
Cross  range  semi  axis  (output) 
Along  track  semi  axis  (output) 
Cross  track  semi  axis  (output) 

None 


semi  axes  by  a  scale  factor. 
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Module  Name: 


Rotate_xy2xpyp 


Calling  Module: 


E_Gate_Semiaxes  3 


Called  Modules: 


None 


Inputs: 


ang_rad 

Z 


Outputs: 


Globals: 


Rotation  angle  (rad) 

Input  two  component  vector 

Rotated  two  component  vector 

None 


Description: 


Rotates  the  two  component  vector  Z  through  angle  ang_rad. 


Equations: 


0  =  ang_rad 


Let  Rot 


_  r  cos(^)  sinf^ 
[-  sin(^)  cos(^ 


Z_p  =  Rot  *  Z 
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Module  Name: 
Calling  Module: 
Called  Modules: 
Inputs: 


Z_m 

Z_P 

ang_rotn 

semi_along 

semi_across 


Test_E_Gate 

Run_Tracker_8 

Rotate_xy2xpyp 


Measurement  state  vector  (x;y;pDop) 
Prediction  state  vector  (x;vx;ax;y;vy;ay) 
Rotation  angle  wrt  radar  cs  (radians) 
Semi  axis  length  along  rotated  x 
Semi  axis  length  along  rotated  y 


Outputs: 

In_E_Gate  Inclusion  test  flag  ( 1  =>  inclusion) 

Globals:  None 

Description: 

Tests  if  measured  point  Z_m  lies  within  an  ellipse  which  is  (1)  oriented  at  an  angle 
“ang_rotn”  with  respect  to  the  radar  coordinate  system,  (2)  centered  on  the  predicted 
position,  Z_p,  and,  (3)  has  semi  axes  of  length  semi_along  and  semi_across.  Returns  a 
value  of  zero  if  test  fails,  and  a  value  of  one  if  test  passes. 


Equations: 

Let  (Xm,  Ym)  and  (Xp,  Yp)  denote  the  measured  predicted  points,  respectively.  Transform 
the  predicted-measurement  position  error  into  the  local  coordinate  system  of  the  ellipse, 
centered  on  the  predicted  point.  The  transformed  residual  errors  (AX',  AY')  are 

[f,;]=Rot. 

where 

Ye  =  semLacross  *  -y/l-(AX7semi_along)^ 

The  condition  for  inclusion  of  (Xm,  YnO  within  the  ellipse  is 

\AY]  <  y,  and 
IAX*]  S  semi_along 


X, -X„ 

Y, -Y. 
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Module  Name: 


T  est_Smile_Gate 


Calling  Module;  Run_Tracker_8 

Called  Modules:  None 

Inputs: 


T_s 

Z_m 

Z-P 

Z_s_last 

x_last_km 

y_last_km 

ag_HI 

semi_rng 

semi_crng 

Outputs: 

In_gate_m 


Scan  period  (sec) 

Measurement  state  vector  (x;y;pDop) 
Prediction  state  vector  (x;vx;ax;y;vy;ay) 
Smoothed  state  vector  (last  update) 

X  position  (last  update) 
y  position  (last  update) 

Maximum  turn  acceleration  (g  units) 
Measurement  ellipse  semi  axis  along  range 
Measurement  ellipse  semi  axis  cross  range 


Inclusion  test  flag  (1  =>  included) 


Globals: 


None 


Description: 

Tests  if  measurement  ellipse  intersects  centripetal  maneuver  “smile”  shaped  gate. 
Equations: 


The  smile  gate  envelope  is  defined  by  two  equations  that  are  parameterized  by  the  time  t, 
ranging  from  the  current  report  time  to  one  scan  period  later.  Let  v  denote  the  constant 
target  speed  and  p  be  the  radius  of  curvature  of  the  turn.  Then 

p  =  vV(9.8  *  ag_HI)  and 
Ang  =  V  (T_s  -  t)/p. 

With  these  the  equations  for  the  maneuver  envelope  in  the  local  track  coordinate  system, 
centered  on  the  last  smoothed  position,  are 

x' =  vt  +  /7sin(ang) 
y.'  =P(l-cos(ang)) 


The  gate  test  is  performed  by: 
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(1)  Transforming  the  boundary  points  (x  s  ,  y  s  )  to  the  radar  coordinate  system 

[x ,  1  r xjast_kml  ^  f cos(e )  -  sin(0 )T x' 

[yj  [y_last_kmj  [sin(e^J  cos(e,ri,)  J[y3_ 

where  0trk  denotes  the  instantaneous  track  angle. 

(2)  Transforming  (Xs,  ys)  to  local  coordinate  system  centered  on  the  measured  point  and 
oriented  at  the  angle  of  the  measured  point  0m- 

'AX^n  rcos0„  sin0„Tx3-x„ 

_ AT/'J  [-  sin0„  cos(0„  )J[y3  -  y„ _ 

(3)  Test  each  maneuver  envelope  point  for  inclusion  within  the  measurement  ellipse 

I  AX 'I  <  semi_rng 

|Ant<Y, 

Y,  =  semi_crng  .y/l-(AX;7semi_mg)^ 
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Module  Name: 


Pred_Shadow_Test 


Calling  Module:  Run_Tracker_8 

Called  Modules: 

get_def_Cov_manu 

predict_shlx_x 

Test_In_Shadow 

Inputs: 

trk_type  Type  of  tracking  filter  (1:3) 

Z_s_last  Smoothed  state  vector  (last  update) 

Cov_s_last_vec  Smoothed  covariance  array  (last  update) 
dT  Time  interval  (sec)  since  last  update 

Outputs: 

Z_p  Prediction  state  vector 

In_Shadow  Shadow  status  flag  (1  =>  in  shadow) 

Globals: 

Cov_man20  Cov_man30 
x_sh_LO  x_sh_HI  y_sh_LX)  y_sh_HI 

Description: 

Tests  if  predicted  position  falls  within  a  shadow  region.  Sets  flag  to  1  if  in  shadow.  [The 
shadow  zones  were  specified  as  a  set  of  rectangular  tiles  by  the  scenario  generator.  The 
array  x_sh_LO,  x_sh_HI,  y_sh_LO,  and,  y_sh_HI  specify  low  and  high  positions  of  each 
time.] 
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Module  Name: 

Calling  Module: 
Called  Modules: 


Test_In_Shadow 

Shadow_Test 

None 


Inputs: 

Z_p 

x_sh_LO 

x_sh_HI 

y_sh_LO 

y_sh_HI 

Outputs: 

In_Shadow 

Globals: 


Prediction  state  vector  (x;vx;ax;y;vy;ay) 

Array  specifying  low  x  position  of  shadow  tiles 
Array  specifying  high  x  position  of  shadow  tiles 
Array  specifying  low  y  position  of  shadow  tiles 
Array  specifying  high  y  position  of  shadow  tiles 


In  shadow  status  flag  (1  =>  in  shadow) 
None 


Description: 

Tests  if  predicted  position  falls  within  any  of  the  shadow  tiles  that  were  specified  by  the 
scenario  generator. 
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Module  Name: 

Calling  Module: 
Called  Modules: 


Discrete_Test 

Run_Tracker_8 

E_Discrete_Test 


Inputs: 

Z_p 

gate_case 

ang_pred 

semi_mg_T 

semi_crng_T 

ang_trk 

semi_trk_T 

semi_ctrk_T 

Outputs: 

In_Discrete 


Prediction  state  vector 
Gate  case  flag  (1:3) 

Angle  (rad)  of  predicted  point 

Semi  axis  length  of  along  range  ellipse 

Semi  axis  length  of  cross  range  ellipse 

Angle  (rad)  of  track 

Semi  axis  length  of  along  track  ellipse 

Semi  axis  length  of  cross  track  ellipse 


Discrete  capture  flag  (1  =>  discrete  present) 


Globals: 


None 


Description: 

Tests  if  discrete  point  falls  within  either  of  two  gates  oriented  along/across  range  or 
along/cross  track. 

Equations: 

(gate_case  =1) 

Each  discrete  point  is  tested  for  inclusion  in  a  single  ellipse,  centered  on  the  predicted 
point,  and  oriented  along  range/cross  range.  The  semi  axes  are  given  by  semi_mg_T  and 
semi_cross_T. 

(gate_case  >1) 

Each  discrete  point  is  tested  for  inclusion  if  either  of  two  ellipses,  both  centered  on  the 
predicted  point.  The  first  ellipse  is  the  same  as  defined  above.  The  second  ellipse  is 
oriented  along  track/cross  track  and  has  semi  axes  given  by  semi_trk_T  and  semi_ctrk_T. 
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Module  Name: 

E_Discrete_Test 

Calling  Module: 

Discrete_test 

Called  Modules: 

Test_E_Gate 

Inputs: 

Z-P 

ang_rotn 

semi_along 

semi_across 

Prediction  state  vector  (x;vx;ax;y;vy;ay) 

Rotation  angle  (rad)  of  ellipse  gate  wrt  radar  cs 

Semi  axis  length  of  gate  along  rotated  x 

Semi  axis  length  of  gate  along  rotated  y 

Outputs: 

In_Discrete 

Discrete  inclusion  flag  (1  =>  discrete  in) 

Globals: 

Ndiscrete  X_D_Vec  Y_D_Vec 

Description: 

Tests  if  any  of  a  set  of  discrete  points  falls  within  an  elliptical  gate:  (1)  centered  on  the 
predicted  point,  (2)  oriented  at  an  angle  “ang_rotn”  wrt  radar  coordination  system,  (3) 
having  semi  axes  lengths  of  semi_along  along  the  rotated  x  axis  and  of  length 
semi_across  along  the  rotated  y  axis. 

The  discrete  locations  were  specified  by  the  scenario  generator  and  given  here  by  the 
arrays  X_D_Vec  and  Y_D_Vec. 
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Module  Name: 


TR_Assoc_Max_T 


Calling  Module:  Run_Tracker_8 

Called  Modules:  cmp_track_age 


Inputs: 

Corr 

Dist 

NATJn 

Nrep_in 

scan 


Binary  correlation  matrix  (trks,  reps) 
Distance  matrix  (trks,  reps) 

Number  of  active  tracks  in  Table_T 
Number  of  reports  from  current  scan 
Current  scan 


Outputs: 

Update 

unAssoc 

unused 

Asgn_Reps 


Binary  array  of  updated  tracks  (1  =>  updated) 
Binary  array  of  unassociated  tracks 
Binary  array  of  unused  reports 
Array  of  reports  assigned  to  each  track 


Globals: 


Table.T 


Description: 

Assigns  unique  reports  to  corresponding  track.  If  a  report  is  common  to  multiple  tracks 
then  it  is  assigned  to  the  oldest  track. 
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Module  Name: 

Calling  Module: 
Called  Modules: 


TR_Assoc_Min_D 

Run_Tracker_8 

cmp_track_age 


Inputs: 

Corr 

Dist 

NATJn 

Nrep_in 

scan 


Binary  correlation  matrix  (trks,  reps) 
Distance  matrix  (trks,  reps) 

Number  of  active  tracks  in  Table_T 
Number  of  reports  from  current  scan 
Current  scan 


Outputs: 

Update 

unAssoc 

unused 

Asgn_Reps 


Binary  array  of  updated  tracks  (1  =>  updated) 
Binary  array  of  unassociated  tracks 
Binary  array  of  unused  reports 
Array  of  reports  assigned  to  each  track 


Globals:  Table_T 


Description: 

Assigns  unique  reports  to  corresponding  track.  If  multiple  reports  are  common  to  a  given 
track  then  the  closest  is  assigned  to  the  track. 
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Module  Name: 


cmp_track_age 
Calling  Module:  TR_Assoc_Max_T 

Called  Modules:  None 

Inputs: 

trk  Track  index  into  Table_T 

scan  Current  scan 

Outputs: 

age  Track  age  since  first  became  firm 

Globals:  Table_T 

Description: 

Computes  age  of  a  track  in  scan  units.  Age  is  defined  as  time  since  track  became  firm. 
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Module  Name: 


get_track_ang 
Calling  Module:  Run_Tracker_8 

Called  Modules;  None 

Inputs; 

7_p  Prediction  state  vector 

Z_s_last  Smoothed  state  vector  (last  update) 

Outputs: 

ang_trk  Angle  of  track  (radians) 

Globals;  None 

Description: 

Computes  track  angle  in  radians  with  respect  to  the  radar  coordinate  system. 
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Module  Name: 


Calling  Module: 
Called  Modules: 
Inputs: 

Z_p 


get_pred_ang 

Run_Tracker_8 

None 

Prediction  state  vector 


Outputs: 

ang_pred  Prediction  of  Angle 

Globals:  None 

Description: 

Computes  angle  of  predicted  point  with  respect  to  the  radar  coordinate  system. 


Module  Name: 

Calling  Module: 
Called  Modules: 


update_error 

Run_Tracker_8 

None 


Inputs: 

N_Hits 

Z_tru 

Z_p 

TID 


Number  of  times  track  captured  a  report 
True  state  vector 
Prediction  state  vector 
Track  ID  number 


Outputs: 

mean_error  Position  error  (predicted  -  true)  averaged  over  N_Hits 

sigma_error  Standard  deviation  of  position  error 

Globals: 

SUMl_ERROR 

SUM2_ERROR 


Description: 

Computes  running  mean  and  standard  deviation  of  position  error  between  predicted  and 
true  values  as  a  function  of  track  ID  number.  Stores  running  first  and  second  moments  in 
global  buffers  SUMl_ERROR  and  SUM2_ERROR. 
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Module  Name: 

Calling  Module: 
Called  Modules: 


Load_P_DET 

Run_Tracker_8 

None 


Inputs: 

Ntk_max 

scan 

N_HITS 

N.EVENTS 


Number  of  tracks  in  P_DET  buffer 
Current  scan  index 
Number  of  captures  for  track 
Number  of  captures  plus  misses  for  track 


Outputs:  None 

Globals:  Table_T_P_DET 

Description: 


Computes  fraction  of  captures  for  specific  track  over  its  evolution.  Stores  results  in 
global  buffer  P_DET. 


69 


Module  Name: 


Load  dBuf_4 


Calling  Module: 
Called  Modules: 
Inputs: 

Ntrk 

dBuf_last_cnt 

Outputs: 

dBuf  cnt 


Run  Tracker_8 


get_Table_Trk_4 


Current  number  of  tracks 

Number  of  dropped  HP  tracks  in  buffer  dBuf  at  last  scan 


Number  of  dropped  HP  tracks  in  buffer  dBuf  currently 


Globals: 


dBuf 

Table_T 

DROPPED  TENT  HRM 
DPLP 


Description: 

Updates  the  high  priority  dropped  track  buffer  each  scan.  The  buffer  dBuf  has  three 
components:  (time,  x  position,  y  position) 
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Module  Name: 


plot_dBuf 

Calling  Module:  Run_Tracker_8 

Called  Modules:  None 

Inputs: 


dBuf_new  Counter  in  dropped  track  buffer 

plot_vec  Min  and  max  s  and  y  values  of  plot  space 


Outputs: 


None 


Globals:  dBuf 


Description: 

Plots  dropped  high  priority  tracks  each  scan. 
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